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Abstract 

A detailed analysis of damping and noise due to a sd-interaction in a thin ferromagnetic film 
sandwiched between two large normal metal layers is carried out. The magnetization is shown to 
obey in general a non-local equation of motion which differs from the the Gilbert equation and is 
extended to the non-adiabatic regime. To lowest order in the exchange interaction and in the limit 
where the Gilbert equation applies, we show that the damping term is enhanced due to interfacial 
effects but it also shows oscillations as a function of the film thickness. The noise calculation is 
however carried out to all orders in the exchange coupling constant. The ellipticity of the precession 
of the magnetization is taken into account. The damping is shown to have a Gilbert form only in 
the adiabatic limit while the relaxation time becomes strongly dependent on the geometry of the 
thin film. It is also shown that the induced noise characteristic of sd-exchange is inherently colored 
in character and depends on the symmetry of the Hamiltonian of the magnetization in the film. 
We show that the sd-noise can be represented in terms of an external stochastic field which is white 
only in the adiabatic regime. The temperature is also renormalized by the spin accumulation in 
the system. For large intra-atomic exchange interactions, the Gilbert-Brown equation is no longer 
valid. 

PACS numbers: 75.10.Jm, 75.30.Gw, 76.60.Es 
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I. INTRODUCTION 



The need for ever higher storage densities and faster retrieval data rates in magnetic 
recording is bringing out new fundamental physical challenges to the industry. In any 
physical device, the main issue is the signal to noise ratio (SNR). For nano-devices, we 
expect a decrease in the signal output and an increase in the noise. Therefore any simple 
scaling-down of the current devices is bound to fail. Hence the need for an understanding 
of the noise source in these devices so that novel solutions to the SNR problem can be 
devised. In magnetic transition-metal based devices, the conduction electrons are one such 
source of noise. The stochastic Landau-Lifshitz-Gilbert equation, which for short we will 
call the Gilbert-Brown equation (GB) 1 ' 2 , has been the main tool in studying the noise at 
the phenomenological level. 3 The GB equation has been quite successful in predicting the 
right form for the damping term for most of the measurements in this area. The noise term, 
which is represented by a stochastic term, is however very qualitative since any meaningful 
account of the noise in a magnetic system is dependent on the microscopic mechanisms that 
give rise to the damping term in the first place. The GB equation has the simple form 



with iieff the effective field, a is a damping constant and the stochastic field h(t) satisfies 



Given that more details about the physics of the devices are now needed to better control 
them, a more microscopic treatment of the noise is in order. This is the subject of this 
paper. However, we will not be able to treat this question in its full generality simply because 
specific details about all systems of interest differ from one to the other. In this study, we 
focus mainly on a thin magnetic film geometry embedded between two normal conductors. 
Such a geometry happens in , e.g., a read- head in a recording device and it is also widely used 
in spin-momentum transfer problems. Damping and noise due to conduction electrons is 
expected to be of importance in these systems. To study this case, we need to compute the 
effective action for magnons in real time for a thin magnetic film embedded in a conducting 
paramagnet. We will derive and solve the corresponding equations of motion, then we will 
discuss the noise spectrum. This program has already been carried out successfully in an 
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earlier paper 4 , where an exactly solvable Caldeira-Leggett-like model has been discussed. In 
that simple case we were able to compute the effective action exactly. However, the model 
of ref. 4 is quite phenomenological and not very predictive, since it involves a very large 
freedom in the choice of the coupling constants. The noise expressions in 4 will however 
surface again in our discussion in the adiabatic limit. 

In more realistic microscopic models, the effective action cannot be computed exactly, 
nevertheless it can be computed in an approximate way. In this paper we will discuss a very 
simple microscopic model, which still contains the basic physics of electrons and magnons in 
thin films and presents quite non-trivial features. This model has been investigated previ- 
ously by many authors, however the relevant references discussing the physics of thin films 
that are directly related to the work presented here are 5-11 . These latter works are pri- 
marily interested in effects similar to the spin-momentum transfer problem of Slonczewski 12 
who studied the influence of a nonzero polarized current on the dynamics of the magneti- 
zation in thin film multilayers. In this paper, we mainly focus on the single film case with 
and without a biased voltage. The case with a polarized current and non-collinear magne- 
tization will be briefly treated numerically in the last section. We will be mainly studying 
the finite size effects of the film on the dynamic of the magnetization but our method will 
allow us also to address questions related to the validity of the GB in the atomistic regime, 
an area of growing interest in recording physics. 13 ' 14 

In the presence of a s-d interaction, the conduction electrons exercise an effective torque 
on the local magnetic moment which can be put in a stochastic form. Therefore, in addition 
to the usual thermal magnetic noise, there will be an additional component due to the 
conduction electrons and one of our tasks is to find out when this contribution can be 
absorbed in the usual Gilbert damping term . The origin of damping in ferromagnets is still 
an open problem. In iron, it is believed that the conduction electrons through the exchange 
interaction are the main channel for the dissipation. 15 In Nickel and Cobalt, the spin-orbit 
coupling is suggested to be the mechanism for the dissipation. 16 ' 17 The calculations of 15 
and 16 are however not totally self-consistent; an adjustable parameter, the relaxation of 
conduction electrons to the lattice, is needed for a meaningful result for the dissipation and 
no treatment of noise has been attempted. Both mechanisms give a Gilbert-form for the 
damping. In this paper we study damping in thin films within the sd-model of ref. 15. 
It has been recently argued that the damping should reflect the geometry of the sample 
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and hence the damping should have a non-Gilbert tensor form. 18 The linear model treated 
in ref. 4 showed that the symmetry of the Hamiltonian has no effect on the damping. 
In this paper, we treat a non-linear interaction between the conduction electrons and the 
magnetization, the sd-exchange, and show that in this case the damping is sensitive to the 
symmetry of the Hamiltonian only for high frequencies. Therefore for the macroscopic 
average magnetization of the sample, the Gilbert damping is correct. Symmetries are 
important only for microscopic magnetization. 

As a result of the recent illuminating work of Simanek 11 , our work will turn out also 
to be intimately related to the spin-pumping theory of Tserkovnyak, Brataas and Bauer 
(TBB) 6 that treated similar questions using scattering theory. Their damping is nicely 
expressed in terms of the mixing conductance, a quantity that needs to be computed by 
ab-initio calculations. In this work we use a very different method which will enable us to 
treat simultaneously atomic magnetic moments and macroscopic magnetic moments simul- 
taneously. Moreover, we will be able to give explicit expressions for the damping and noise 
at all frequencies and include finite size-effects of the film within the sd-model. Realistic 
systems can also be treated by this method but will require numerical computations. Hence 
our results will be of interest to those interested in atomic simulations of magnetic systems, 
an area which is starting to become important for magnetic recording. 14 

The experimental work of Covington et al. on spin momentum transfer is another mo- 
tivation for our work. 19 This latter work showed that in a biased spin valve with currents 
below the critical current, i.e. current needed to switch the thin layer, the resistance shows 
large 1/f-type noise in the MHz-GHz regime. Our system is similar to a spin valve except 
that we do not have a reference layer. This will enable us to examine the contribution of the 
sd-exchange to the line- width in the N/F/N structure and the spin momentum transfer. A 
full micromagnetic treatment is also given that includes the effect of non-spin flipping events 
on the spin-momentum transfer. Our conclusions will be helpful to the interpretation of 
the experiment and the micromagnetic calculations. We will show an example where the 
noise has its origins in the non-uniformity of the in-plane magnetization. 

The paper is organized as follows. In sect. II, we set up our notation and the Hamilto- 
nian used in our calculations. We use a non-isotropic Hamiltonian that takes into account 
the ellipticity of the magnetization which is typical in thin magnetic films or local magnetic 
moments. Using the real-time formalism for our model, we first compute the free propa- 
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gators of the theory. Then, we derive the effective action of the system by integrating out 
the electron degrees of freedom. We derive a stochastic equation for the magnetic moment 
that is different from the GB equation. These equations are non-local in space and in time. 
These equations will be especially needed in atomic-type simulations of magnetic systems 
where the local effective field is large compared to macroscopic fields. In sect. Ill, we 
discuss the limits under which we can recover the GB equation for this model in a macro- 
scopic system N/F/N. We show that interfaces enhance the damping and the fluctuations 
of the magnetization. The bulk damping is assumed to be due to conduction electrons 
interacting with the lattice. In sect. IV, we show how to calculate the noise spectrum 
of the magnetization and discuss its dependencies on the anisotropy, initial conditions and 
on the spin accumulation for the N/F/N case. Moreover, we show that in the adiabatic 
limit this sd-exchange interaction is equivalent to a stochastic external field with a Gaussian 
white noise distribution and effective temperature that reflects the geometry of the system. 
For high frequencies, we show that the GB equation is no longer valid and that the damping 
reflects the symmetry of the Hamiltonian which does not appear in the linear or adiabatic 
regime. In sec. V, we discuss a geometry similar to that of ref. 19 using a macroscopic spin 
transfer model. Based on the quantum calculations in previous sections and the micromag- 
netic calculations, we suggest that the noise in 19 is due to thermally assisted transitions 
between two non-uniform states of the magnetization. Finally in the appendix, we give 
various expressions needed in the calculation of the correlation functions of the magnetic 
moment and discuss the dependence of the damping on the symmetry of the Hamiltonian 
of the magnetic moment. 



II. HAMILTONIAN FORMULATION AND THE ONE-LOOP EFFECTIVE AC- 
TION 

Let us consider a thin film of magnetic material interacting with a large external magnetic 
field of order one Tesla or more H = (0, 0, H) , constant in time, uniform in space and 
directed along the z direction which is in-plane. Let us assume the film has linear dimensions 
DxL y xL z , i.e. it has a rectangular section in the yz plane (the plane parallel to the magnetic 
field), with area L y L z and thickness D in the x direction, with D « L y and D « L z , as 
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shown in figure 1. The magnetic film will be later assumed to be sandwiched between two 
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FIG. 1: The geometry of the thin film described in the text. 



large normal conductors. We are interested in studying the effects of the exchange magnetic 
field due to conduction electrons on the average magnetization of the thin film as well as 
the local atomic moments. Ultimately, we will be interested in studying the case D is much 
smaller than the lateral dimensions of the film, i.e. the thin film limit. The finite size effects 
related to the breaking of translation invariance in the x direction will be of primary interest 
to us since in this case the effect of the conduction electrons on the average magnetization 
is strongest. A path integral formulation proves to be very useful in problems of this sort. 

A ferromagnetic material is roughly a system of quasi-free electrons of spin s (the 4s 
electrons) interacting with bound electrons of spin S (the 3d electrons) via an Heisenberg 
Hamiltonian 



where J is the interaction constant, of the order of at least 0.1 eV in the macroscopic case 
and about 10.0 eV in the microscopic case, a = (<7i, a 2 , a 3 ) is the vector with components the 
Pauli matrices and ip is the 2-component electron field. The 3d bound electrons are mostly 
aligned with the external magnetic field and they are the source for the magnetization S of 
the material. The Hamiltonian of the conduction electrons is then given by 
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FIG. 2: The local magnetization vector. 



where the potential V(x) is a spin independent confining potential of the structure. For a 
large external field H, the 5' 2 -component of the magnetization can be taken to be a constant 
and hence its interaction can be absorbed in the diagonal part of the energy. 

The simplest possible effective Hamiltonian we can construct for the magnetization is 



-h • s + l -c xx s x {*) 2 + \c yy sy^f 



(3) 



where c is the cell index with volume v = V/N. In micromagnetic simulations, v is of the 
order of 10 nm 3 while in atomistic calculations it is of the order of the 0.4 nm 3 . We have 
neglected higher powers of the S* and components. This can be justified once we notice 
that in a magnetic material the spins are mostly oriented in the direction of the magnetic 
field, i.e. the spin vector has locally x and y components which are small with respect to 
the z component: 

s c = (s:,sy,s z c ), \s*\,\sv\<:s* (4) 

In other words, in this paper we will consider the case in which the spin vector has a small 
angle with the z axis (see Figure 2). Micromagnetic calculations show that an angle of 20 
degrees can still be considered small . The dynamics of large angles and in particular the 
possibility of magnetization switching is quite interesting too, but it cannot be addressed 
within the approximations used here. 20 

The advantage of the small angle approximation, i.e. taking S* time-independent, is the 
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simplification of the commutation relations, S x S = iS, 21 ' 22 which become 

[S* c , S y c ] = tS z c S c c • (5) 

This implies that S y c is canonically conjugate to S*. We neglect any exchange stiffness 
between the cells. Next, we consider the case of an infinite wavelength spin wave with 
wave-vector k — 0. This corresponds to an homogeneous spin field. The Heisenberg 
equation of motion can be derived from the Hamiltonian and the commutation relation and 
they assume the form 

S x CyySy , Sy Cxx^x • (6) 

The coefficients c xx and c yy are related to the anisotropies of the medium and we are espe- 
cially interested in the strongly anisotropic case, c xx ^> c yy . The Heisenberg equations of 
motion can be trivially solved in the form 

= e M(t-to) 

\Sy(t) J 

where M. is the matrix 

M=(_° C "\. (7) 
More explicitly, we have 

S x (t) = coscjqO - t ) S x (to) + J^- sm u (t - f ) S y (t ), 

V Cxx 

Fc 

S y (t) = ~\ — sino;o(t - *o) S x (t ) + coscu (t - t ) S y (t ) , 
V °yy 

where o;o = ^Jc xx c yy is the frequency of the elliptic precession of the magnetization. In the 
rest of the paper, we calculate the effect of the conduction electrons on this solution in both 
low and high frequency limits. Since the interaction with the non-dynamical S z component 
has been already accounted for in the Hamiltonian of the conduction electrons, only the a\ 
and o"2 terms enter in the interaction with the dynamical spin field. In momentum space 
this interaction term reads 





n 



int 
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J 

qk 



K-k^i«q^k + a[-^2a q Sl] . (8) 



The interaction Hamiltonian Timt commutes with the total Hamiltonian. This means 
that only two interactions are possible: a) a spin up electron +1/2 emits a +1 magnon 
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and becomes a spin down — 1/2 electron; b) a spin down —1/2 electron absorbs a spin +1 
magnon and becomes a spin up +1/2 electron. These interactions are represented in figure 
3. 



a) 



b) 



FIG. 3: Conservation of angular momentum in the z direction 

Next, we calculate the effective action of the theory to first order in J 2 . This is done 
by eliminating the electronic degrees of freedom after which, we get a stochastic equation 
for the magnetization valid in the non-adiabatic regime. 

First we compute the free propagators which are obtained by inverting the diagonal 
differential operators id t — for the electron field and the differential operator 



5S 

1 " 6Si k (t)SSi(t>) 



■xx &t 

(I Cyy 



5 P (t-t') (9) 



for the spin field. Here S k and S are the diagonal parts of the action for the electrons and 
the magnetization, respectively. Therefore, one has to solve the differential equations for 
the corresponding propagators Qp and Dp: 

(id t -£ k )g P (t,t',k) = 5 P (t-t') (10) 

and 

V- 1 V P (t,lf,k)=8p(t-lf) (11) 

where the delta function is defined on the closed-time path (CTP) P and the boundary 
conditions are such that the spin field is periodic on the path P while the electron field is 
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anti-periodic. 23 Notice that eq. 10 is a 2 x 2 matrix in the spin space whereas eq. 11 is a 
2x2 matrix in the complex plane. The first equation is diagonal and can easily be solved 
with solution 

g p (t - t') = -ie-^-^Qpit - t') + Ae'^-^ 

where the matrix A is the integration constant corresponding to a generic solution of the 
homogeneous equation. The boundary condition Gp(t ,t') = —e^^Qp(t — i(3,t') fixes A = 
ifk- Therefore 

Q P (t, t', k) = -i (0 P (t - f ) - f k ) e -^-*') . (12) 

The second equation 11 seems more complicated, but actually can be simplified by mul- 
tiplying both sides by ia 2 and using the identity ia 2 V~ l = d t — M. , where M. satisfies 
M T = —a 2 M.a 2 . A simple computation gives 

V P (t, t') = e M ^ ia 2 [Q P (t - t') + n(tM T )} e MT{t '^ o) , (13) 

where the integration constant matrix n{iM) is fixed by the boundary conditions Vp{t , t') = 
V P (t -if3,t') to be 

n{iM T ) = 1 . (14) 

exp(i(3M T ) - 1 

Notice that the time evolution is trivial and an explicit computation gives 

( sin uj ( t -t') cosuj (t-t r ) \ r 
° V } ° V } \[&p(t-t')+n(tM T )} (15) 

— cos u (t — t') — ^ sin u (t — t') I 

Since M. has eigenvalues +iu and — iuo, the eigenvalues of n(iA4 T ) are n(— cj ) and n(u ) 
respectively, i.e., they are regular Bose-Einstein distributions. A little algebra allows to 
derive the free spectrum as 



<: {S*(t),Si(t)} :>= 



Cyy sinh Quip _^ 
wo 1— cosh/3u>o 

j Cxx sinh §wp 

wq 1— cosh f3wo 



In the isotropic case, things are simpler since M. is proportional to ia 2 and therefore 
it commutes with ia 2 . We also observe that in the CTP formalism, each of the matrix 
equations 12 and 13, corresponds to an additional 2x2 matrix of equations in the Schwinger 
formalism, depending on the position of t and t' in the path Pi directed along the positive 
evolution in time and P 2 directed in the opposite direction. 24 In particular 
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g P (t,t') = G>(t,t') 



1 _ / k ) e -**(t-f) , teP 2 , t'eP 1 



Q P (t, t') = G<(t, t') = ife-^W , teP 1 ,t'eP 2 



(16) 
(17) 



It is easy to check that the previous propagators are consistent with the operator's expres- 
sions 

G>(t, f) = -i< a k (f)4(f') >c G<(t, t')=i< a£(t>k(f) > c (18) 
= -i < SL v {t)S>tf) >c = -i < Sj(O5ik(0 > c • (19) 

Next, we include the effect of the conduction electrons on the magnetization. Using 
the functional formulation, 24 it is clear how to extract the effect of the electrons on the 
local magnetization field: it is enough to integrate out the fermionic degrees of freedom 
and to compute the effective action for the magnetization. Since the original action is 
quadratic in ip and the functional integral can be performed. If we use a coherent state 
representation for the magnetization and the conduction electrons, the generating functional 
for the problem has the form 25 

Z[J] = J dz x dzl f dd.dQ, exp [- Zl zl - Ci,iCi,i] -Ci \pU ] ( J 2 ) U (JO KCi) (20) 

where z's represent the d-electrons degrees of freedom while £'s represent the conduction 
electrons. J\ and J 2 are the usual virtual sources. The propagator inside is given by 

(||) = J dzdz* J d(dCexp [z\z (t - t(3) + Cr >fcj C fcj (to - ifi) (21) 

rto—i/3 

i I ds(iz*d s z-H J (z*,z)) 



to 



+i 



where 



(k, P ) 



[id s -e 1 (k))5(k-p) 
iz* (k - p) 



iz(k-p) 



(id s -e 2 (k))S(k-p) 



(22) 



the variable s represents the time along the CTP P. In the above we have set the total 
number of magnetic cells N — 1 since we the same discussion applies to a local atomic 
moment. Next we integrate out the conduction electrons degrees of freedom. Using 
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J d(*d(exp i J 



exp [-Tr In G" 1 ] , 



(23) 



dsCG-\ 

and expanding the logarithmic term gives exchange terms of all orders in the coupling con- 
stant J. In the following, we keep only quadratic terms. The conduction electron propagator 
satisfies the equation 

'(id s - e t (k)) 5 tJ 5 kM - V l3 (k - k')] Q%' (s - s) = -8 ip (s - s) 5 (k - k') (24) 
with V, the interaction term is given by 

z (k - k') 
z* (k - k') 

This equation is solved by iteration, assuming that higher order terms are small. The large 
Zeeman-type term has been included in the e term. The propagator Q is then given in 
terms of the propagator which is a solution of the following equation 



V( k - k ')-| 



(25) 



s, s 



-Sij [s — s ) 5 (k — k j 



— is \ s — s 



(26) 



(id s - e> (*)) Q\T 
If we define the functions G < and G > to be 

G< ( S ' S ) = _ ^ £ ) exp 
G> (s, s') = i (1 - /(e)) exp -is (s - s') (27) 

then formally, we have 

(s, s') = G< (s, s') 0p (V - s) + G> (s, s') 0p (s - s') , (28) 

where 6p (s — s') is the step function along the closed-time path P. Hence, the presence of 
conduction electrons in the thin film will give rise to an effective action on the d-electrons. 
To order J 2 , the logarithmic term becomes 



Tr InQ = \ £ J 2 f dsdsg^ k (s, s) (V) ^ («\ *) V^ k (s) (29) 

fc,fc' 

= J 2 ^ / dsds'z (s' , k — k' j 2* ^s, k' — kj exp — i (^Ei (k) — e 2 ( s — s ) 

fc,fc' c 

x [6 (a - 5) / (e 2 (k')) (1 - / ( £l (k))) + 9 (V - s) f (e 1 (k)) (l - / (e 2 (k')))~ 
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which can be written in a more compact form as follows 



Tr In g = J 2 [ dtdt exp -% [e x (k) - e 2 (t - t' 

k) % (t-t, k, k' ) Zj (t', k - k') 



(30) 



k,k 

x z* (t, k' 



In matrix form, 



a n (t-t',k,k') = e(*-t')/(e 2 (k'))(i-/( £l (k))) 

+0 (t'-t)/( £l (k))(i-/( ea (k'))), 

2t 12 (W,k,k') = /( ei (k))(l-/(e 2 (k'))), 
2l 21 (t-t',k,k') = /(e 2 (k'))(l-/( ei (k))), 
2t 22 (W,k,k') = e(*'-*)/(e 2 (k'))(l-/( ei (k))) 

+0 (f-f')/( ei (k))(l-/(e 2 (k'] 



(31) 
(32) 
(33) 

(34) 

where / (e CT ) is the Fermi-Dirac distribution for spin-up (er = 1) and spin down (er = 2), 
respectively. The effective action of the magnetization is now given by 



S e " = J dt [iz*d tZl -H J (zlzj] 

- J dt [iz* 2 d t z 2 - H J (z*, z 2 )] + iTr In Q 
where z — S + — S x ± iS y . Next, we make a change of variables 



(35) 



z 2 ) \l -1/2 ) \ 3 
The new quadratic form is then 



= T 



3 

3 



(36) 



(37) 



This form first appeared in Schwinger's paper but is often cited under the name of the 



13 



Keldysh form 26 , 



U 
Un 

Wi 2 

W21 

W22 



fi-i',k,k' 
(t- t',k, k' 
(t-i',k, k' 
fi-i',k,k' 
ft-f',k,k' 



T T 21T, 


0(t'-t) [/( £l (k))-/(e 2 (k'] 
-0 (t-t) [/( ei (k))-/(e 2 (k' 
\f ( £l (k)) (l - / (52 (k'))) + i/ (.2 (k')) (1 - / ( £l (k))) . 



(38) 
(39) 
(40) 
(41) 
(42) 



It is easy to see from these definitions of the kernels Uij that 

Wi 2 (k,k',t) = -U 21 (k,k',-t), 
U 22 (k,k',0=W22 (k,k',-t). 



(43) 
(44) 



We can now introduce a new field 77 into the theory. This will mimic a random Gaussian field 
in the semi-classical equations. In the generating functional Z [J] , we replace the quadratic 
term in 3 and 3* by linear terms, 



cxp 



^ fcfc' 



exp 



-2 ( £l (k) -e 2 (k')) (*-*')] 3 (V,k-k') 



/ dr]*dr]exp - ^ / rft / rft V (*> k ) 25-1 (*> k ! *7 (*'» k ' ) 

^ \_ kk' J J 



(45) 



where the kernel D is given by 

J 2 



D t 



,k,t , k ^ 



-i ( ei (k)-e 2 (k')) (t-t') 
x (/( £l (k))-/(e 2 (k')))coth 



(46) 



'(5 ( ei (k)-e 2 (k')) 



In the bulk, this noise kernel will in general depend on the relaxation of the conduction 
electrons due to phonons which will appear in the exponential term on the right hand side. 
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After this integral transformation, the effective action becomes 

-^(3*3 + 3*3)-^(3V + 33)} 

+ [ dt [ dt ' ex p H ( £i ( k ) - £2 ( k )) ^ - 
\ ' kk' J J 

x {3* (t, k' - k) W 12 (t-t, k, k') 3 (*', k - k') 

+3 (t, k' - k) U21 (t-t', k, k') 3 (V, k - k') } 

+ J dt{v* (t,k)i(t,k)+i*(t,k)ri(t,k)} + F[ri*,rj\. 



(47) 



The equations of motion are obtained by minimizing the action with respect to the variables 
3 and 3*, 

5 (S°") 

= 48 



Si (*,p) 



r)*=r)=Q 



and 



V (t, x) = Vx (t, x) + ir] y (t, x) (49) 

We find that the transverse components of the local magnetization satisfy the following 
equations of motion 



dS x (t, p) 
dt 



(n-K)S y (t,p)-rj y (t,p) 



+ (£) S / d * ,/C ( tj k; k + p ) cos ( k ) " £ 2 (k + p)) (t - t')] S- x (t, -p) 
+ (0 S / d * ,/C (*' k ' k + p) sin ( k ) - £ 2 (k + p)) (f - f')] S y (t, -p) , 



(50) 



and 

ffi/ (*» P) 



-(n + K) s y (t, p) + rj x (t, p) 



[ dt }C (*> k 5 k + p) sin ( k ) - £2 (k + p)) - *')] 5 X (Y, -p) 
+ (£) S / d * ,/C (*' k; k + p ) cos (k) - £ 2 (k + p)) (t - t')] ^ (t, 



-P 

(51) 
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where 

K (t, k; t, k + p) = 6 (t - t) [f ( ei (k)) - / (e 2 (k + p))] (52) 
and the correlation function for the random field due to the conduction electrons is given by 

V (t, p) rj (7, P ) ) = y (/ (ei (P)) - / (e 2 (p'))) coth ^ 
x exp -i (ex (p) - e 2 (p 



2 l £ MP 



ei IP, 



(53) 



The correlation function of the x-component is then given by 
ite (ifc (f, -k) ifc (f\ k') ) = (^) (/(ei(k))-/(e 2 (k')))ooth 



cos 



(W) 

(54) 



where o; fc)fe / = e\ (k) — e 2 (k')- This is one of the important results in this work. The kernel 
terms account for the dissipation and a shift in the frequency due to the interaction with 
the conduction electrons. These equations of motion differ from the usual Gilbert form since 
the dissipation term is not a derivative and is non-local. These equations generalize those 
derived by Mills for a Stoner particle at zero temperature and in the adiabatic limit. 9 Hence 
the memory terms will be very important for local moments. In the next section we study 
the limit under which the GB equation is recovered in a thin film embedded between two 
large reservoirs at equilibrium. 



III. FINITE-SIZE EFFECTS IN THE MEMORYLESS LIMIT 

In a series of very illuminating papers, Simanek was able to show how the ideas of 
Tserkovniak, Brataas and Bauer can be understood in the familiar linear response approach 
which avoids the use of the scattering method. 6 ' 8,10 Similar calculations were carried out by 
Mills 9 using a dynamic RKKY approach which generalizes the earlier results obtained by 
Berger 5 . It is well known that the spin-pumping theory and the Berger-Mills theory both 
give interfacial additional damping due to spin currents (not charge currents). In the former 
theory, this damping vanishes when there is no exchange splitting between the spin-up and 
spin-down electrons while it does not in the Berger-Mills theory. The spin pumping theory 
however seems to be very successful in interpreting the recent experiments by Mizukami et 
al. 27 . Hence in this section, we use the equations derived in the previous section to further 
understand this particular discrepancy between the various methods. Our results happen 
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to be similar to those derived by Simanek using the spin-pumping theory. We believe 
however that our approach is more direct and transparent besides it is self-contained. This 
equivalence has very important consequences on the understanding of the physical origin of 
the spin momentum torque in finite films. Moreover, our theory is easy to extend to finite 
temperature and can deal with transient conditions as we show in the next section. 

First, let's set-up the geometry of the problem and calculate the damping in the limit 
when there is no memory in the magnetic system, i.e., the average magnetization is much 
slower than the conduction electrons. This is the adiabatic limit. The geometry we adopt 
(fig. 4 is the same as the one adopted by TBB 6 . The two reservoirs on each side of the 
thin film will act as a sink for the spin leaked through the interfaces. The reservoirs are 
maintained at the same chemical potential fi in this section and hence there is no net flow of 
charge from left to right. Our theory can be also adapted to the case of non-equal chemical 
potentials which is briefly addressed in the following section. 



M 








z 











-D/2 D/2 



FIG. 4: A thin film confined between two large reservoirs with the same chemical potentials /x = //. 
The case of different /x's is treated in sect. IV. 



For the rest of this section we use Mills 9 notation since he was able to derive a more 
general form for the Gilbert equation in the adiabatic limit. We will show below how our 
theory reduces to his in this limit. 

Within the linear response approach, the Gilbert equation for the magnetization M is 

dM 



dt 



[Mx(H + (H e// (t)»] 



G\l\ 



M x 



dM' 

~dt 



(55) 
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where G is the Gilbert constant and the effective field H e ff is due to the action of the 
conduction electrons with the magnetization through a sd-type interaction 

H'ff(t\ - j2Vc (a dM * + a dM y\ 



H*ffm j2Vc (a dM y A dMx \ 



AG = W^2^\^ M s- (57) 



Hence the conduction electrons enhance the Gilbert damping term by 

J 2 V C 
2Nhfi 2 

It was moreover argued by Mills that the constant Ai, which renormalizes the precessional 
frequency, is not zero in general as was assumed by Simanek and Heinrich 8 . Below we show 
that our analysis naturally gives an explicit expression for this term and that it vanishes 
within the approximations employed here. In higher orders in J 2 , the contribution of this 
term is non-zero but small as we show in the next section. 

From the equations of motion for the x,y-components, Eqs. 50-51, we easily see that it 
is the term that has the sine-dependence that gives rise to damping, 

72 poo 

^G x(y) (t, p) = - / dt'J2 [/(ei(k) - / (e 2 (k + p))] (58) 

x sin (ei(k) - e 2 (k + p)) t'] S y{x) (t-t, -p) . 

The cosine-dependent term will be shown below to be the origin of the term Ai discussed 
by Mills. If we use the slow (adiabatic) precession approximation for the spin variables, 

S x(y) (t-t, -p) = S x(y) (t, - P ) - /W*»-P) + ... (59) 
Inserting this back in Eq. 58, we get 

AG x(y) (t, p) = AG^ y) (t, p) + AG<J tf) (t, p) , (60) 

with 

72 poo 

(*> P) = (*' E / dt> U( £ ^ - f (e 2 (k + p))] 

xsin (ei(k) -e 2 (k + p))*'] (61) 
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and 



AG d v) (t,p) = 



J 2 dSy (x) (t, 



dt 



roc 

^[/( £ i(k)-/( £2 (k+ P ))] 

L, JO 



xt'sin 



(ei(k) -e 2 (k + p))*' 



(62) 



The term AG^X contributes to the precessional frequency while AG^ gives rise to Gilbert 
type damping. Hence it should be related to the imaginary part of a susceptibility term. In 
fact we can write 



AG« } (t,p)=-< 
where the 'susceptibility' function x is 



■ dx (fi,p) 



^(a) ~p) 
dt 



X (t, p) = -y 6 (f) £ [/( £l (k)) - /(e 2 (k + p))] sin [f ( £l (k) - 5 2 (k + p))] , 

k 

and its Fourier transform 

/oo 
X (*, P) e^di 
-oo 

= -xE^( k )-/^(k+ P ))] 

k 

1 



(63) 



(64) 



x 



+ 



(65) 



_£i(k) - e 2 (k + p) + Q + it] £i(k) - e 2 (k + p) - Q - ir\_ 
where 7] is a small positive real number. From this expression, we see that to allow for a 
finite relaxation time r s of the conduction electrons, we can just replace Q by Q + i/r s . The 
constants A 2 introduced by Mills can be obtained from x through the following expression 



dtt 



= A 2r (p) +iA 2 (p). 



(66) 



An easy calculation shows that A 2r vanishes within this approximation and A 2 , which is still 
a function of the momentum, is given by 



A 2 (p) = -J 2 ^[/( £l (k))-/(£ 2 (k + p))]( £l (k)-e 2 (k + p)) 

k 

7] 



(67) 



(( £l (k)-£ 2 (k + P )) 2 + 7 7 2)^ 

As shown by Heinrich, Fraitova and Kambersky 15 in an infinite medium, the conduction 
electrons can't dissipate energy unless the electron-hole pairs have a finite-lifetime by trans- 
ferring energy to the lattice. This can be done by taking a temperature-dependent finite 
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rj. In this case the damping term A 2 (p) will not vanish. This spin-flipping mechanism 
is believed to be the source of the damping in iron and permalloy 28 . Another source for 
damping can be geometrical in origin. As shown by Mills a breakdown of wave vector 
conservation due of the finite size of the film can give rise to dissipation. In the language 
of reservoirs, we rephrase this by saying that quenching the states of the magnetic film to a 
countable number while leaving those of the electronic states denumerable is equivalent to a 
Caldeira-Leggett (CL) model with a fermionic bath which is known to give rise to quantum 
dissipation. 4 ' 29 Our geometry is then a typical example of this model and should show a 
dissipative behavior as a function of the thickness D, i.e., A 2 (p) — > as D — > oo, even in 
the case of very slow relaxation times for the conduction electrons. We take the thin film 
to have finite thickness D in the x-direction, and we take the transverse components, S x 
and S y to dependent only on the x-coordinate normal to the plane with pinned or unpinned 
boundary conditions. In the continuum approximation, the magnetization components in 
this symmetric configuration take the form 

Si(t,r)= Sr(i) cos (^) , i = x,y (68) 

ji=0,1,2,... 

The Fourier transform is given by 



dx Si(t,x)e tpx = Si(t,p) (69) 
therefore, we have Si (t, —p) = S* (t,p) and 



S l (t,p) = (2nf5{p ll )J2s:(t) 



sin ((p + f ) D/2) sin ((p - f ) D/2) 



(70) 



where we have set p x = p, — oo < p < oo. 

In the following we assume that the splitting 5 in electronic energy bands due to the 
sd-interaction is smaller than the Fermi energy and k^T « \i which is the case at room 
temperature, then it is enough for our purposes (because of the finite size) to use the following 
approximation for the Fermi-Dirac functions for the conduction electrons 

[/(e 2 (k + p))-/(ei(k)] = |^( £2 (k + p)- £l (k)) 

Bf > (*■»+!), (71) 



where Ef = (e F + £ f) /2- This approximation is not necessary and will not change our 
conclusions but it helps keep the algebra at minimum, otherwise the Lindhard function 
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will appear explicitly in our expressions and will considerably add to the complexity of the 
calculations. If a non-zero voltage difference is applied across the thin film, then Eq. 71 
has to be modified to take account of the spin accumulation effect A u = /i T — ^ due to the 
normal-ferromagnetic interface 30 . This term which can be positive or negative depending 
on the direction of the polarized current will hence contribute to the damping. We will say 
more about this case when we study the noise and the corresponding fluctuation dissipation 
theorem in the next section. Using the approximation Eq. 71, the damping term becomes 
(for p ^ 0) 



A 2 (p) = 




8 — V\ 



5 + p 



in J ~ s 



(s+ P 



P 



2e, 



m 



tan 



5 + p 



2e, 



m 



(72) 



where the Fermi energy ep is that of the spin up electron in the ferromagnet. Care is 
needed to get the corresponding expression for p = 0. The damping is therefore momentum- 
dependent as we should expect in a finite film. This expression can be, e.g., useful for 
studies of spin-wave resonance in thin films. In the remaining, we confine ourselves to the 
volume mode since it is usually the mode measured by FMR. In this case, the damping 
simply becomes 



A 2 ~ 



4 



dp 



sin 



(P?) 



A 2 (p). 



(73) 



'o V V 

As it is clear from this expression, damping is directly related to the breakdown of momentum 
conservation in the direction normal to the film. This damping clearly vanishes when the 
size of the film becomes infinite, i.e., D — > oo. This expression for the damping was 
found by Fourier transforming back the equations of motion , eqs. 50-51, to real-space and 
enforce the condition that the magnetization vanishes for \x\ > D/2. This is the procedure 
that we followed to allow us to capture the finite size effects of the magnetic film on the 
damping in this section and on the noise in the next one. Using the fact that the saturation 
magnetization is defined by M s = /Iq/v, the excess Gilbert damping is therefore given by 



G s = ^A 2 . 

v 



We write this in terms of small dimensionless parameters 5 r 



— and r r 



(74) 

]rjsF With 
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typical values 0.4 and 0.001, respectively in transition metals. The damping is now given by 

27r 2 hD~ 



Gs = *rM^F(T r ,5 r ,D), (75) 



with 



F{T r ,o r ,D) = — / ax— sin -kfDx 



4 ; ^[( 5r _ ;r ) 2 + T 2] 

<5 r + x tan -1 [(5 r — x) /r r ] 



x [(5 r + x) 2 + r r 2 ] xr r 



In the limit of vanishing bulk damping, r r — > 0, and large thickness D, the damping is given 
by 

_ fi 2 k 2 F f 0.555 r \ 
Gs -2^h{-^~)- (77) 

To get this result, we used the approximation 

f°° , s ^ °- 55 

for large x. Hence this damping is due to finite-size effects. At smaller thicknesses, the 
damping shows some oscillatory behavior as shown in figs 5,6,7 for small bulk damping. 
The size 10 8 sec -1 ) and the oscillatory behavior of the damping near the interface are 
similar to what Mills 9 and Simanek 11 found and is non-existent in a large sample as is clear 
from the figure for large D. The dimensionless Gilbert damping for iron would be therefore 

a = — w 0.01, (79) 

which has the right order of magnitude as measured in ref. 27. Figure 8 is perhaps the case 
that applies to transition metals. In this case the oscillations are almost nonexistent and is 
consistent with the recent numerical calculations of Zwierzycki et al. 31 which are based on 
the circuit theory approach. 32 The temperature dependence is weak in all the results since 
we have assumed that ksT << ep. For large thicknesses, the damping is therefore still 
dependent on the exchange coupling and this dependence is linear. This is in contrast to 
Berger's 5 result where his interfacial damping is independent of J. 

Finally, we show that the term Ai vanishes to lowest order. To find an explicit expression 
for the Ax-term, we proceed along similar lines as we did for A 2 . We make a slow time 
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FIG. 5: The damping constant (in sec 1 ) as a function of the thickness of the ferromagnetic film 
(in nm). k F = 10 8 cm" 1 , S r = 0.2, rj r = 0.001 
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FIG. 6: The damping constant (in sec 1 ) as a function of the thickness of the ferromagnetic film 
(in nm). k F = 10 8 cm" 1 , S r = 0.1, m = 0.001. 

approximation for the cosine-term in eq. 50 and we define a 'susceptibility' function x' p) 
similar to x, 

X ' (n, p) = -iJ 2 (f (ei (k)) - / (e 2 (k + p))) (80) 

k 

1 1 

X 

Ei (k) - e 2 (k + p) + Q, + it] ei (k) - £ 2 (k + p) - ^_ 
As in the Ai case, A 2 is proportional to the imaginary part of dx' (fi = 0, p) /dVt which is 
easily seen to vanish. At higher orders in J 2 it gives a nonzero contribution to the frequency. 
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FIG. 7: The damping constant (in sec 1 ) as a function of the thickness of the ferromagnetic film 
(in nm). k F = 10 8 cm" 1 , S r = 0.2, rj r = 0.001 
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FIG. 8: The damping constant (in sec 1 ) as a function of the thickness of the ferromagnetic film 
(in nm). k F = 10 8 cm" 1 , S r = 0.6, rj r = 0.001 

This is the subject of next section. 

IV. THE MAGNETIC NOISE SPECTRUM 

In this section, we calculate the various correlation functions of the magnetization vector 
by including higher order corrections in J in the exchange field. We will deal with both 
the low frequency limit and the high frequency regime. The former is applicable to the 
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case of large magnetization while the latter is important for local atomic moments. We 
will show that the damping for local moments depends on the symmetry of the Hamiltonian 
and will not be of the Gilbert-form. In the adiabatic limit which applies to the average 
magnetization of the film, the GB equation is recovered. 

This is a direct extension of the calculations presented in the previous sections. In the 
bulk it was shown by Heinrich, Fraitova and Kambersky 15 that the sd-exchange gives zero 
damping unless dissipation to the phonons is included. This was done in a non-self-consistent 
way by putting a relaxation time r s by hand in the electron propagator. However, based 
on a simple analogy with Migdal's theory 33 on electrons and phonons, we should expect a 
nonzero damping for spin waves with wave numbers p such that < k F — kp < p < k F + kp 
if ksT « e F (fig. 9). Hence the volume mode, i.e, p = mode, won't be expected to show 
any dissipation in the bulk. Sine we are partly interested in long-wavelength excitations in 
thin films, we will concentrate on the consequences of interfacial effects which based on the 
calculations above should provide a new source for dissipation. We follow closely Schwinger's 
original work on the harmonic oscillator interacting with a bath of harmonic oscillators. 24 
Some of the computations will be deferred to the appendix in the hope not to distract the 
reader from the end results. Our treatment is self-contained. Similar calculations have 
been carried out in Ref. 4 which dealt with a hypothetical physical model for dissipation 
in the bulk. However in the adiabatic limit, we will recover the results in 4. This shows 
that the macroscopic magnetization as opposed to the local moments is insensitive to the 
dynamics of the environment. In this section, the bulk spin relaxation time r s = 0. We 
go beyond the mean field approximation for the exchange field. 34 In other words, we seek 
corrections to the propagators of the theory by including self-energy corrections. This will 
allow us to go beyond the approximations made in the previous sections and calculate the 
damping due to inelastic scattering of the conduction electrons off the magnons as they 
cross the interfaces. We do not use the slow-time approximation in this section, but the 
damping will be shown to have the Gilbert form for frequencies much smaller than the 
electronic precessional frequency. Moreover, we find a very interesting result that relates 
the GB equation for spin-momentum transfer to the model treated here. We show that the 
Langevin dynamic treatment of Li and Zhang 35 for the noise can be only justified for the 
adiabatic limit, and around the FMR frequency in thin films. As in the previous section, 
the magnetization precesses around an effective in-plane field directed along the z-axis. We 
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FIG. 9: A spin-down/spin-up excitation in a bulk ferromagnet with smallest and largest wave- 
vectors. 

start by coupling the transverse magnetization S = (S x , S y ) to external sources Ji and J 2 
along the positive and negative-time oriented paths, respectively. The modified Hamiltonian 
is now given by 

Wp = W 1 -J 1 .Si-(W 2 -J a -S 2 ). (81) 

Next we introduce the new variables J = | (Ji + J 2 ) and Q = Ji — J 2 . Similarly, we define 
an average S = | (Si + S 2 ) and a difference D = Si — S 2 representing the fluctuations in the 
magnetization. The generating functional Z [Ji, J 2 ] is defined in such way that its Taylor 
expansion at Ji = J 2 = gives the correlation functions of the magnetization, 

(82) 



(83) 



(^pT 1 exp 












Texp 


-i [ f dtH[J 2 ] 


) 








J to 





Then we see from the definitions that 

51ogZ[Ji,J 2 ] 



and 



* <Si(*)> 



Ji=J 2 =0 



i{Si(t)Si(t)) 



(84) 



The latter symmetric average is the one usually associated with the noise in the magneti- 
zation vector. We will seek a general expression for this quantity that takes into account 
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initial conditions, i.e., the reservoirs and the film are assumed initially to be separately in 
equilibrium before they are put in contact with each other at time to- Hence our method is 
capable of handling transient behavior in our system. To calculate the symmetric correlation 

function, we need first to find the equation of motion of the average value (S) . For this, 
we need the effective action in the presence of the external fields which is easily found from 
Eq. 47 



^ e// [3*,3,3*,3] = £ r 



dt 



~y (k) da (k) - \? (k) $3 (k) 



+ ^3 (k) d t f (k) + 1 3 (k) d t y (k) - in (3* (k) i (k) + 3 * (k) 3 (k)) 
-iKtf (k)3*(-k)+ 3 (k)3(-k))] 

+ J' 2 J2 [ dtdt'e-^W^ [3* (t, p) U 12 (t - t') 3 (f, -p) 
k,k' •* 

+r (t, p) u 2 i (t - f) 3 (t, - P ) + f (t, P ) u 22 (t, f) 1 (?, - P )] 

+ i J2 [ dt (Q* c 3 + Q c 3*) + i J2 \ dt ( J *d + J c3*) , (85) 

where we have p = k' — k and e 12 (k, k') = £x (k) — e 2 (k'). The remaining variables are 
defined as follows: 

3 = ^(5, + <5 W ) , (86) 

3 = ^= (A. + <I>y) , (87) 

Q c = (Q x + iQ tf ) , (88) 

Jc = 4= ( J - + * J J/) • ( 89 ) 

v2 

Varying the action with respect to its variables, we get the respective equations of motion 
for the magnetization and the fluctuations. A solution of these equations will require a 
careful treatment of the boundary conditions. To do this, we found it easier to take a 
related representation of the coherent states. This representation is equivalent to the usual 
harmonic oscillator representation. 

We define two harmonic oscillator-type operators a and a + such that 
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Hence, the energy operator for the spin system becomes 

H = fajjQ (a + a + . (91) 

with the frequency given by 

CJ = Sy/CxxCyy J ( S = 1 ) . (92) 

Therefore coherent states in this representation are defined by 

a\Z)=Z\Z). (93) 

To change between representations (3,3) — ► (Z,z) ( with careful handling of the boundary 
conditions) in the path integral representations we need to make the change of variables 

3 -> \ (ciZ + c 2 Z*) , 

3 -> ^ (ci* + c 2 2;*) , (94) 
with the coefficients C\ and c 2 given by 



Cl = 



C2 = — v^—' (95 » 



It will be seen that in the high frequency regime, the damping depends separately on c\ and 
c 2 . The path integral representation on the closed-time path requires doubling of variables 
as has been done above. The following boundary conditions are also needed: 

Z 1 (t f ) = Z 2 it f ) , (96) 
Z 1 (t ) = exp [-Pu ] Z 2 (t ) , (97) 
Z* (t ) = exp [+pu ] Z* (t ) . (98) 
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It is important to observe that the boundary conditions are non-hermitian in the coherent- 
state formulation. This is a direct result of the application of the Kubo-Martin-Schwinger 
(KMS) condition. This implies that Z, the average sum of Z x and Z 2 , and their respective 
fluctuations z satisfy the conditions 

Z (to) = ~ ooth (^) z (to) , (99) 

and 

z(t f ) = 0. (100) 
The free action in the new representation therefore takes the form 

<s<°> [z; Z, z\ z\ = \ {Z f Z} (t f ) + Z}Z f (t f )) + \ (Z,Z* (t ) + Z*Z (to)) 
+ \{Z iZ * (to) + Z*z(t )) 

+ ij\t{^Z*z-Z*z) + ±(z*Z-z*Z 
-lu (z*Z + zZ*) + Z*A + ZA* + z*Y + zY*} , (101) 
where the external sources A, E are defined through the relation 

3 1 ■ Si - J 2 • S 2 = Z*A + ZA* + z*Y + zY*. (102) 

The equations of motion are obtained from the full action iS e ^[Z*,Z,z*,z] = 
^g(o) [z*,Z,z*,z] + iS B [Z*,Z,z*,z], where the last term is due to the interaction of the 
conduction electrons with the magnetization. For the fluctuations z and z* , we obtain 



and 



(d t + icu ) z(t,p) / dt'G 1 (t - t', p) z (if, -p) 
-Jlfdt'G 2 (t-t',p)z* (t',-p) 

= tA (t,p), 



(dt - iw ) z* (t, p) +^ / dt'Gz (t - f, p) z* (If, -p) 
+£fdt'G 2 (t-t',p)z* (t',-p) 

= -iA*(f,p). 



(103) 



(104) 
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To solve these equations and take account of the finite size of the film, we introduce the 
following definition for the 'averaged' Green functions, 

— i r r r D/2 r D/2 r d 3 v 

G a {t-if) = ^J d\\ J All J dxe^ J dye*™ J j^G a (t - t>, p) . (105) 

where A is the area of the surface of the film. This definition is useful for solving for the 
zero mode only. As we did in the previous section, this is obtained by Fourier transforming 
back the equations of motion to the real space representation where the condition of finite 
thickness can be easily implemented. These Green functions (a = 1, .., 5) are given explicitly 
in the appendix. It should be noted that the Green functions G Q with a = 1,2,3 are directly 
related to the appearance of the dissipative term in the magnetization. They involve terms 
similar to those that appeared in Eq. 71, but they are also the Green functions that appear 
in the equations of motion of the fluctuations z and z*. Hence their dissipative nature is 
very clear in this formalism since they introduce irreversibility in the dynamics of Z and Z*. 

The corresponding Green's function needed for the solution of the fluctuations are there- 
fore given by 

(d t + iu ) Sl (t-t')-^J dfG x it - t") 01 (*" -t')=5(t- f) , (106) 

01 (t - t') = 0; t> t' (107) 

and 

(d t + iu ) fli (t -t') + — / drGi (t - t) ; (r - f) = 5 (t - t') , (108) 

q[ (t -t') = 0; t< t' (109) 

Similarly, we get two more equations for Z(t) and Z*(t) that involve two more Green 
functions. We only write the respective solutions below. Therefore, the solutions of the 
fluctuations and the average sum take the form 



z(t)=i J & (t - f) A [t^j dt 

~ 1 (£) / dTdt ' dt "^ (* - *') ^ ( f - *") & (*" - r ) a * - ( n °) 
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and 

Z(t)=i J dt'g[ (t - t') E (f) +i^J drdr'dt'^ (t - r) G A (r - *') 01 (t' - r) A (r') 
_ i 3 ! J drdr'dUx (t ~ t) G 5 (t - t') q\ (t' - r) A* (r 1 ) 

-i^J drdT f dt'Q , 1 (t - t) G 2 (r - t') & (t' - r) S* (/) . (Ill) 

Since z(tf) = 0, we require that 0i (t — t') = for t > t'. 

First, we recalculate the initial correlations to show that we have the correct boundary 
conditions. The initial value for Z follows from the solution for Z(t) after setting the coupling 
constant J = 

Z (t ) = ~ coth (/^) i J dt% (t - A (*') . (112) 

Then the derivative of the x-component of the magnetization with respect of the external 
sources is 

i&) cotb ("t) 101 (i " + 9; (t " 01 (113) 

In the free theory, the propagator jji is simply given by 

Sl ( t - t') = -0 it' - t) e — «(*"*'). (114) 

Hence 

-—coth f/3 — Icoswni 
2w„ V2 / 



i 5 A J s (t) 2 cj 

= l({S x ,S x (t)}), (115) 

which is the desired relation that was derived in Sect. II by a different method. To get 
this solution, it was crucial that we apply the correct boundary conditions on Z and z, Eqs. 
96-100. 

For the coupled case, the initial condition for the Z(t) equation is 
Z (to) = ~\ coth (p^p) | -i dt' 9l (t - t') A (f ) 

- i (£) jH drdt'dt" Ql (t - t') G 2 (t - t) Ql (t" - t) A* (r) . (116) 
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Hence the general solution for Z(t) is 



Z (t) =Z (t)+i / dr dl (t - t) E (r) 
•/to 

/ I {] [ drdr'dt'g, (t - t 1 ) G a (t - r) 01 (if - r') A (V) 

2 <>oo 

/ drdr'dt'g, (t - r) G 5 (r - £ (f - r') A* (r') 

2 />oo 



(117) 



2 

""I 



— i 



J 



drdrdt'^ (t - r) G 2 (r - *') (*' - r') £* (r') 



(118) 



to 



where Zq (t) is a particular solution of the nonhomogeneous problem 



-iwo(t-to) _|_ ( ^ ) p -iuot f di'pi-^ot' 



z (t) = z Q (t) e -^-^> + j - | < 



to 



'G 2 (f-f )e UJ °y t ~ to )Z (t ) 



I dt " 

Similarly, we get the solution for Z*(t) 

POD 

z* (t) = z* (t ) e ^ o(t - to) - % \ dt'oi (t - t') s* (0 

•/to 

- * (0 / drdr'dt'Q* (t - r) G 2 (r - t') g[ (f - r') £ (r') 
/ drdr'dt'g* (t - r) G 5 (r - t') 01 (*' - r') A (r') 



</M"m (/ t )e-^°( t "- t0 )z (t ) 

(119) 



7 t 
2 coo 



+ i (£) jT drdr'^Vi (i - r) G 4 (r - fit (*' - r') A* (r') 



J dt'G^t-t^e^'-^Z^to) 



-J dt"G 2 (*'-r) e -K'''-<°)zo(*o) 

with the initial state given by 

Z (t ) = ~i coth y di'fl! (t - A (f) . 



(120) 



(121) 



Now, it is easy to calculate the components of the magnetization from the above results. 
We just need to differentiate the average magnetization with respect to (j x ,j y ), 
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V2 




(A + A*) , 



V2 [W (A A *\ 



(122) 
(123) 



to find the symmetric correlation functions for the magnetization. 

A straightforward calculation gives the differential of the x-component with respect to 

Jx 



% 5j x (tf) 2 uj 



|_I coth [ e -o(*-*o) * {tQ _ f) + e -iMt-to) Qi {tQ _ ^] 



e -^ (t-to) / rfr / dr / e ^o(r-t ) 



to 



- e ^ o( ^' o) J dr J dT'e-^ o{T - to) [g 3 (t - r) e ^ o(r '-' (,) g* (i„ - t') 
-G 2 (r - r') e-^'-*°) 01 (t -0]] 

y rfrrfr' [q[ (t - r) G 4 (r - r') 01 (r' - 

- flf (f - r) G 5 (r - r') 01 (r' - f ) - g', (f - r) G 5 (r - r') fl ; (r' - f ) 
+ f (t-T)G 4 (T-r')gt(r'-t')]]} 

= \({s x (t)XV)})- (124) 

This is a general result for the noise due to spin-flip scattering between magnons and 
conduction electrons that is useful for all frequencies. It is clear from this expression, that 
the correlation functions depend explicitly on the initial state of the system. However, 
there is a term which is invariant under time translation. This part of the correlation term 
survives at times much later than the initial conditions. Before we turn to the calculation 
of this term which is the term usually measured in FMR-type experiments we give the 
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expression for the correlation function of S x and S y components of the magnetization 

i r i 



5jy (t) 2 



-- coth [e-^(*-*») fll (t - - e**<*-*> fl ; (f„ - 0] 

gf (f - r) G 5 (r - r') 01 (r' - f ) + * (t - r) G 5 (r - /) £ (r' - if) 



-tf (f - r) G 4 (r - r') fl J (r' - f ) 

2 



3 iwo(T— to) 



x y di" [Gi (r - t") e --°(*"-*°) 01 (t - t') +G 2 (t- t) e^'-^Ql (t - t') 

+ 1 Q 2 e -o( t - t o) CQth ^ dTe -^o(r-to) 

x I dt" G 3 (r - t) e -o(*"-*o) s * (to - f ) + y *"G 2 (r - f") e -*-o(*"-*o) fll (t(J _ t >) 

(125) 

An explicit expression for these correlation functions is not needed for what follows, but 
we will write its limit in the adiabatic limit which is the case of most current interest for 
the average magnetization of the film. Before we do that, we would like to make few more 
comments about these general expressions for the correlation functions. These expressions 
are beyond the usual fluctuation-dissipation relation and hence they can be adapted to 
truly 'non-equilibrium' situations. As an example we mention time-dependent pulse-field 
excitations of the magnetization, switching by a magnetic field. In this latter case, the 
z-axis is the local equilibrium axis. They also give us an idea on the noise behavior in local 
magnetic moments and the corresponding damping. 

To calculate the noise spectrum in the x-component of the magnetization we need the 
Fourier component of the function 

C xx it - t>) = (S x (t) S x (0>l t » t0)t '» to , (126) 

since here we will not address the transient regime, which is treated elsewhere 20 . Its Fourier 
transform is easily found to be 

C -H = ^(£) Re {si M G 4 (u) di (u) - 0i (w) G 5 (u) flj (-co) } . (127) 
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The damping for large thicknesses D (and nonzero exchange splitting) acquires a simple 
asymptotic expression, 

\ 



7T (C\ + 4) ( Sin (l ( k F ~ 4) D) COS U (fcj - 4) D) 

*> — 1 + 2 / y +48 ^ y 7 + ... 

(4 + 4) s ^ (4 - 4) d (4 - 4) D * 



(128) 



This damping differs from the one found in the previous sections in two different aspects. 
The first is that for large thicknesses, the damping is weakly dependent on the sd-exchange 
energy 5 r as opposed to being linear in 5 r . This result is now much similar to Mills and 
Berger. The second important difference is that the relaxation time depends explicitly on 
the symmetry of the original Hamiltonian of the magnetization. This result is however 
similar to what we found in ref. 4. For circular precession, we simply have cf + c 2 , = 1 and 
hence any dependence on the form of the precession is lost. The damping is still however 
a scalar of the Gilbert form for uj « 8/H and does not appear to require a tensor form as 
suggested in ref. 18. The high frequency regime, which is applicable to atomic moments, 
is however more interesting in this respect. The relaxation time is not a simple function 
of the ellipticity and hence the damping is not of the Gilbert form; the relaxation time is 
an algebraic function of the frequency (see appendix). The damping in the film still shows 
oscillations as a function of the thickness of the film and it attains larger values than in the 
previous calculations to first order in J 2 (fig. 10-11). 




FIG. 10: The damping as a function of the thickness for S r = 0.2 and c xx = lOOc^. The last term 
is a measure of anisotropy and is on the high side for realistic permalloy films. 
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FIG. 11: The damping as a function of the thickness for S r = 0.6 and c xx = 100c yy (same as fig. 
10). 

The expression for the noise, Eq. 127, shows a very interesting property that only the 
Green functions G4 and G§ show up as multiples of the coupling constant J 2 in the nu- 
merator. The dissipative Green functions do not appear in this form. This observation is 
important when we try to include the effect of spin accumulation on the noise and under- 
stand why the effective temperature concept appears in the stochastic formulation of Li and 
Zhang 35 for close to equilibrium. It is also important to observe that in the half-metallic 
limit, i.e., J — > 00, we have 

Hence, the exchange energy splitting drops out completely from the damping expression. 
This result is similar to that derived by Bazaliy, Jones and Zhang 36 in a half-infinite fer- 
romagnet in contact with a normal metal. In the adiabatic limit, u « J,the expression 
for the noise can be simplified considerably. At high temperature (i.e., k B T >> w and 
T «T C the critical temperature), it becomes 

_ 2ak B T ( uj 2 + Cyy + {auf \ 

l c i + c 2) \ [{oj - loq) + [aw) \ {{lu + Lu ) + {aw) \ J 
where the damping a is thickness dependent and vanishes when D — > 00 in the film. The 
shift in frequency is less than one percent and is neglected in the macroscopic case. A 
similar result derived in the isotropic case has been recently communicated to us by 37. A 
comparison of this expression with the corresponding expression in ref. 3 reveals a very 
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interesting result; the dynamic of the d-electrons in the presence of the sd-exchange inter- 
action is simply reproducible by the GB equation with white noise and a is replaced by 
a G ( c i + c l) where «g is the Gilbert damping. This is an important simple result which 
shows the reasons behind the successes of the GB equation. 3 At the same time this latter 
result points to the limits of applicability of the GB equation in atomic simulations. 14 The 
recent work of Safonov and Bertram 38 suggests changing the damping form of the GB equa- 
tion. Our calculation clearly shows that for the average classical magnetization this is not 
needed. It is only for high frequencies that the effect of the symmetry of the Hamiltonian 
on the damping becomes appreciable (see appendix). For low frequencies, the dependence 
of the relaxation time on the ellipticity factor, cf + c^, agrees with Kambersky and Patton 39 . 
In this work, we are also able to give explicit expressions for the damping term and the 
corresponding noise in thin films. In the next section, we will use the GB equation with a 
damping a = 0.03 to study the continuum case. We will also include polarization due to 
another magnetic layer with relative angle far from zero or 180°. Next we treat the effect 
of spin accumulation on the noise in the low frequency regime and close to the FMR fre- 
quency. As we stated after Eq. 71, the spin accumulation term will have to be included in 
the self-energy term. The only Green functions where this has to be taken into account are 
those that appear as a result of the fluctuations z and z*. These Green functions will affect 
only terms where the parameter a and the frequency u> appear in the Gilbert form in the 
G\ Green function only. The overall a-term that shows up in front of the temperature will 
therefore remain unaffected by the spin accumulation term. To show this result, requires 
a more careful evaluation of the Gi(u). The k-integrals over k and k — p are carried out 
separately. In the limit of small damping and close to the FMR frequency, we can define 
an effective temperature in the GB equation 




(131) 



and a renormalized interfacial damping 

a " = a { l + ^p)- (132) 

This will allow us to write the corresponding noise in a form that is remarkably similar 
to that derived recently by Li and Zhang 35 where they considered only bulk damping in 
their problem. However there is a very important difference. The difference in chemical 
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potential is expected to be proportional to the current 5,7 but not linear to the damping as in 
35 which is based on the Slonczweski picture. Therefore this expression should be able to 
differentiate between the sd mechanism and the Slonczewski mechanism. There are already 
indirect indication from ref. 40 that the critical current is not linear with damping which is 
the result derived in 35 and includes only bulk damping. In CPP structures therefore we 
should not expect the spin momentum transfer term to change the character of noise in these 
systems as long as we stay below critical currents. There are however other effects in this 
geometry not addressed by this calculation. One of them is the field from the current and 
the other is magnon-phonon interactions which my may affect the average magnetization 
configuration in these films. We discuss these two effects classically in the following section. 



V. DISCUSSION AND CONCLUSION 



In this final section, we would like to make some comments and discuss our results in 
light of recent experiments on noise in spin valves, such as that in ref. 19, which show 1/f- 
type noise in CPP structure . In CPP spin valves with biased fields, it was observed that 
excessive low frequency noise is generated with current. It was argued that spin momentum 
transfer between the magnetic layers is responsible for this noise. In this section, we treat 
a case where the spin momentum transfer is not the root cause of noise in these geometries 
but it will just affect the amplitudes of the noise. The noise in our case will be inherently 
due to at least three contributing factors: the biasing of the spin valve, the field from the 
current and the thermal fluctuations in the system. We reach these conclusions based on our 
calculations carried out in previous sections and on simulations based on the GB equation 
with a spin torque as it was suggested in ref. 12, 

^ = 7 S x (u eff + a G ^ + h) + PIS x (S x S p ) , (133) 
with the stochastic field h(t) satisfying 

{h^hjit')) = 2a G k B T5 lj 5{t - t'). (134) 

The parameter f3 is a geometrical factor and I is the current. The pinned layer with 
magnetization S p is not dynamical and hence the finite size effects discussed above can 
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be included by choosing a large ota- The effect of spin accumulation is not taken into 
account properly in the stochastic field; we are simply assuming the effective «g to be a 
constant and current independent. This equation is solved numerically with a G = 0.02, 
M s = 1400 emu/cc and thickness d = 3 nm. 41 ' 42 The white noise approximation as discussed 
above is valid for frequencies around the FMR frequency which is of the order of 10 GHz. 
There are two theories of spin momentum transfer: one is microscopic and based on the sd- 
exchange model while the second is macroscopic and is based on a simple balance equation 
for the spin currents. The calculations presented in previous sections are closer to the first 
approach rather than to the second one. It is believed that the first approach is dominant 
only at very thin films of 10 nm or less. 43 Therefore based on the results derived here, it 
appears that any 1/f-type noise measured in CPP spin valves should be attributed to non 
sd-type of interactions. We show below that the spin torque does not appear to generate 
1/f-type noise and it is the non- homogeneities in the magnetic configurations that are mainly 
responsible for the noise. The CPP structure we study is shown in fig. 12 and is similar to 
that in 19. The current is flowing from the pinned to the free layer which are separated by 
a normal conducting layer. The single particle simulations of the magnetization show no 
interesting behavior and are noiseless and this is consistent with the results from the previous 
section. The field from the current is taken into account in the calculations and is needed to 
observe 1/f-type noise at finite temperature. The magnetization is also biased in the y-axis 
and x-axis with a 300 Oe field and a —90 Oe field, respectively. The demagnetization field 
of the sample is also taken into account. The magnetic material is chosen to be that of a 
permalloy. Figures 13 and 14 clearly show that the effect of the spin torque only slightly 
increases the already present noise in the system for the particular parameters shown in 
the figure. It does not give rise to the low frequency in this example. A closer study 
of this example shows that it is the combination of the biasing, the field from the current 
and the temperature that are the source of the noise. Hence in this example, neither the 
sd-type model nor the macroscopic model can explain the origin of the 1/f-type noise. A 
calculation that does not include in great detail the configuration of the magnetization is 
therefore highly unlikely to capture the source of the noise in these structures. Figures 16 
and 17 show the two possible metastable states that are responsible for the 1/f-type noise 
in this device. 
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FIG. 12: The CPP spin valve: The thick layer with magnetization S p is pinned along the x-axis. 
The magnetization Sf is free to move. The current I is perpendicular to the interfaces. Spin 
momentum is transferred from the pinned layer to the layer by polarizing the current with the 
fixed layer. 
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FIG. 13: The xx-component of the noise spectrum for the average x— component of the magneti- 
zation with the spin torque T sp between the magnetic layers set to zero in the GB equation. 

In summary, we have studied in some detail the noise and the damping problem in thin 
magnetic films embedded between two normal conductors. We have mainly focused on the 
interaction between the conduction electrons and the d-electrons as the main mechanism 
for damping. Our results also apply to the microscopic case at high frequencies. In this 
model, it has been shown that the damping in thin layers oscillates as a function of the 
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FIG. 14: Same as in Fig. 13 but with the spin momentum transfer torque included in the GB 
equation. 
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FIG. 15: A real-time trace of the average magnetization component S x (only a typical interval of 
time is shown). The magnetization appears to oscillate largely in this component for this magnetic 
configuration where the total magnetization is mainly along the y-axis that is perpendicular to the 
polarization axis. 

thickness of the film. If higher orders in the exchange coupling constant are taken into 
account, the relaxation time becomes dependent on the ellipticity of the precession of the 
magnetization. Only at high frequencies, the damping is no longer of the Gilbert form 
and becomes explicitly dependent on the ellipticity. This result is important for atomic 
simulations. The noise associated with the interfacial damping has also been calculated. It 
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FIG. 16: Configuration of the magnetization in state 1. The horizontal arrows are along S p . The 
angle between the external bias field and S p is close to 90 degrees. 







/. 


/. 


/ - 


/. 


/. 


/. 


/ . 












/. 


/. 


/: 


/ 


/ ' 


/ ^ 


t 


/ 


/ 


/. 


/ . 


/. 


/. 


/ 


/ 


i 


I 


/ 


/ 


/ 


/ 


/. 


/. 


/ 


/ 


t 


t 


/. 


/ 


/. 


/. 


/: 


/. 


/ . 


/ ; 


I 


/. 


/. 


/. 


/. 


/. 


/. 


/; 




/ 


t 


/. 


/. 




/ 


/. 


/. 


/. 








/. 












/. 


/. 


/ 




/. 


y. 


y 


y 


y 




y. 


y 


/. 





FIG. 17: Configuration of the magnetization in state 2 (same as in fig. 16). 

was shown that it does not give rise to a large 1/f-type noise. The spectral density curve gives 
the usual peak at the natural frequency of the system at low damping. In the adiabatic 
limit, we have found that for small damping and close to equilibrium, the spectrum is 
well represented by a white noise source term and an effective temperature that measures 
deviations from small spin accumulations. Therefore a large negative current can give rise 
to a large renormalized damping and suppression of the noise amplitude around the FMR 
peak while positive current decreases the effective damping and increases the noise around 
the FMR peak. Finally, we have also shown that even within a macroscopic calculation the 
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transfer of spin momentum does not give rise to large low frequency noise, rather in our 
example it was the field from the current, the thermal fluctuations and a particular biasing 
that simultaneously give rise to large noise at frequencies below the FMR frequency. 
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APPENDIX A: APPENDIX 

Here we give the definitions of the different functions that appear in the main text and 
briefly discuss the high frequency regime of the damping term which becomes dependent on 
the symmetry of the Hamiltonian in a non-trivial way. The following function are derived 
by integrating out the conduction electrons degrees of freedom from the magnetization- 
conduction electron equations of motion. The functions are 



G 1 (t-t',p) = J2{cpi2 (t-t')e 



-iei 2 (k+p)(t-t') 



(Al) 



k 




G 2 (t-t',p)=J2{ 



C1C2U12 (t - t') e 



•iei 2 (k+p)(t-f) 



(A2) 



k 




and 



G 3 (t - t', p) = J2 { c ^i2 (* " e-^^+PW-* 1 ) 



(A3) 



k 
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are due to the correlation terms between the two branches of the CTP path of the path inte- 
gral. The following functions are symmetric in time and will be the origin of the correlation 
functions of the random field: 

G 4 (t - t', p) = J2 { c i U ™ (k + P) e-^+pX'-O (A4) 

k 

+ cjU 22 (k-p) e --( k -P)(*-*')j, 
G 5 (t - t', p) = J2 {^U 22 (k + p) e -^ + m-t') ( A5) 

k 

+ Cl c 2 U 22 (k-p) e «»(k-p)(t-f)|. 

All these functions are not independent. For example, we have 

G 3 (t-t',-p) = -G 1 (t-t',p). (A6) 

In most of our calculations, the functions U 22 are approximated by the following expres- 
sions: 

U 22 (k + p) = l -5 (e (k) - e F ) (^-P + coth + a) , (A7) 

and 

W 22 (k - P ) = ±* (e (k) - e F ) + a) coth + a) . (A8) 

These approximations are not necessary but makes the algebra less involved. 
The function G 1 has the exact explicit expression: 

G 1 (t - t', p) = J2 {" c ? @ (* " [A ( k ) - h (k + p)] e- fe » (k+p)( "') (A9) 

k 

+c\Q (f - t) [fx (k) - j 2 (k + p)] e i£ -( k -P)(*-*') } . 

Similar expressions for G 2 and G3 can be easily deduced from that of G\. 
All these functions are needed for the calculation of the propagators for the fields Z,Z* \z 
and z*, which are to first order in J 2 given, respectively by 
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and 



01 (t -t') = Q{t- t') e^ ^ (AlO) 

+ (€\ r dt " r dre (t"-t) e **(«-«") 



xG, {t"-r)e{t'-r) e^H) 



+ ..., 



'dt"y rfre (Y'-t) e -** (*-'") 

xG 3 (t"-r)e(/-r) e--'H') + ..., 

; (f - t ') = e(t- t') e"^'('-'') (A12) 

"(0 s ' J dt " J dT<d (*-^") e-^°(*-*") 
xG, (t"-r)0(r-t') e -^H') + ..., 

g ; (t - f') = e (t - t') e ^°(*-*') (A13) 



+ ( ^ / / drQ [t-i 



x G 3 (7' - r) O (r - *') e 



'"'II ( ' ; 



The various Green functions needed for the calculation of the noise spectrum are given 
here in Fourier space (h — 1) . The free propagators associated with the fields Z(t) and z(t) 
are 



flf } M = --, (A14) 

UJ — UJq + If] 

C M = (A15) 

— cj — «?y 



In the presence of the conduction electrons they become (A = J/2), 
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and 



si M 



si M 



0i o) M 



1 - A 2 Gi ( W ) gS 0) (w) 



0i (O) M 



(A16) 



(A17) 



We will also need their conjugate form, which is again different from the complex conjugate. 
The free propagator is 

(A18) 



/*(0) , . 

01 M 



and the dressed one is 



0!* M 



u + uj + tr] 

'*(0) / X 

0i M 



(A19) 



1-A 2 G 3 Hfli (0) H' 

The kernels G Q , (a = 1,4,5), are averaged over all spin wave modes. For Gi, it is defined 
by 

2 



d (w) = D 



/• d 3 k 


r + °° dp 


sin (p§) 


/ (27T) . 


l-oc 2^ 


Pf 



(A20) 



x < c 



•9/ /fcxP 



+ c: 



2 df_ f-k x p 
dei 



+ 5 



5 



5 + ^ - uj + 



2?7 



m 



5 - + w - ™ 

with similar definitions for the other kernels. The Fourier transforms are approximately 
given by 

2 

(A21) 



r°°d P 


sin (py) 


Jo 27T 





x <j cf ( 2/c F + — log 



+ c^\2k F + log 



p 
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5 + u + hFP 
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G 4 M 
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sin (p§) 


/o ^P 


P# 



(A22) 



x <^6 A; F 



m (a; — 5) 



p 



+ c^e ( k F 



m{uj + 5) 
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and 



— m A vc x c 2 / u\ 

G5 M = — 2 — ucoth {^2) 



rdp 


sin (pf ) 


to n P 


P? 



(A23) 



x < 6 A; 



m (a; + 5) 



7; 



fc 



m (a; — 5) 



P 



In all of the above expressions, k F = F 2 F ■ It should be also noted that 

G\ (u) = G 1 (-uj). 

Finally we give the term that is directly responsible for the relaxation term, 

2 



(A24) 



TimDv 
/■ceGi = u 



dx 1 sin a; 

x \ x 



clQ 



2ki 
mD 



-x 



(A25) 



where c\ = (y / c^+ y 7 ^) / (4cj ) and c\ = (v^r + v^ra) / (4^o) • For the noise in 
the average magnetization, the range of frequencies we are interested in are usually very 
low compared to the exchange splitting energy. In this case the relaxation time which is 
proportional to TZeGi will depend only on the overall ellipticity factor c\ + c\. The cases of 
higher frequencies are of interest only in the atomistic limit which will also give a similar 
expression as in A25 for the relaxation time. In this case, the relaxation as easily seen 
will depend separately on c xx and c yy and is no longer linear in uj. Memory effects in the 
systems become important and the GB equation is no longer valid. In the linear model of 
ref. 4, we did not have any higher order dependence on frequency and the damping was of 
the Gilbert form. Hence representing a bath by harmonic oscillators is only useful when we 
are interested in the low energy limit. 
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